Detailed analysis of quantum phase transitions 
within the u{2) algebra 



L.Fortunato^Q and L.Sartori^ 

1- Dipartimento di Fisica "G.Galilei" and INFN, via Marzolo 8, 1-35131 Padova, 
Italy. 

2- CISAS G. Colombo, via Venezia 15, 35131, Padua, Italy 



Abstract. We analyze in detail the quantum phase transitions that arise in 
models based on the u(2) algebraic description for bosonic systems with two types 
of scalar bosons. First we discuss the quantum phase transition that occurs in 
hamiltonians that admix the two dynamical symmetry chains u{2) Z) u{l) and 
u{2) D so(2) by diagonalizing the problem exactly in the u{l) basis. Then we 
apply the coherent state formalism to determine the energy functional. Finally we 
show that a quantum phase transition of a different nature, but displaying similar 
characteristics, may arise also within a single chain just by including higher order 
terms in the hamiltonian. 
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1. Introduction 

Quantum phase transitions are a subject of current interest in various branches of 
physics that range from molecular to nuclear and to hadronic physics [H [21 [3l Hj. 
Often one describes a certain set of measurable properties of a quantum system with 
phases that have a genuine quantum origin and, within some simplified algebraic 
model, tries to associate these regimes with the manifestation of a given (dynamical) 
symmetry. These symmetries are connected with Lie algebras and, by exploiting their 
mathematical properties, one is usually able to make detailed predictions for important 
observables and classify large sets of data using recurring patterns and schemes that 
proceed from the formal algebraic structures. These methods have been heavily 
employed in nuclear and molecular physics, for example in the Interacting Boson 
Model (IBM) and Vibron model[5l[6], that make use of the u(6) and w(4) algebras 
respectively. However, due to the high dimension of the algebras involved, it would be 
easier to discuss the simplest possible case to highlight that several important features, 
discussed in connection with the IBM or other algebraic models, are extremely common 
(and very likely universal) as they appear even in the paradigmatic case of u(2). This 
algebra is often found in connection with the schematic Lipkin model, that is often 
taken as an example of solvable many-body system. Two subalgebras may occur within 
the u(2) algebra: one, u{l) is usually associated with a vibrational spectrum and with 
a 'spherical' phase, while the other, o(2), is associated with a rotational spectrum and 
with a 'deformed' phase. In the following we will study the ii(l)-o(2) phase transition 
(the only possible within this scheme), by solving a transitional hamiltonian in the 
u{2) D it(l) basis. We will give energy spectra along the transition and study the 
phase transition in terms of the mixing of components of the original eigenvectors. 
The analysis of the composition of the eigenvectors reveals a remarkable persistence 
of the vibrational character during most of the phase transition up to the critical 
point. 

In parallel to the analytic calculations that made within the Lie algebraic scheme, 
coherent states have been used to connect the abstract algebraic formalism to a 
geometric description by means of a mapping of second quantization operators in 
terms of differential operators. This procedure allows to transform the algebraic 
hamiltonian into a Schrodinger-like differential equation that can be used to calculate a 
potential energy surface (the expectation value of the hamiltonian in the ground state) . 
The position of the minima of this functional determines in turn several important 
properties of the initial system and furnishes a way to study and interpret the phases 
and the possible phase transitions as well as their dependence on certain parameters. 
The coherent state has been heavily employed in the IBM, in the Vibron model and 
in several other models. This method is crucial to define the critical point and the 
character of the phase transition. 

Moreover, we will use both the exact diagonalization and the coherent state 
formalism to show how a quantum phase transition between two different behaviours 
can be obtained within a single subalgebra chain, just by allowing a combination of 
linear and quadratic Casimir operators of u(l). 

The u{2) algebra has been investigated in the thesis of O.van Roosmalen [7], 
in several books about algebraic models [6] and in a recent paper by Cejnar and 
Iachello[2^, that treats the u(n — l)-o(n) quantum phase transition in general. Our 
paper provides a detailed discussion of the simplest specific case of u{2). 

In Ref. [8] a scalar two-level boson model is set up to simulate the consistent Q 
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formalism (CQF) for the Interacting Boson Model getting rid of the complications of 
the quadrupole boson operators. This allows the model to be solved for a much 
larger number of bosons, that is a severe limitation in the usual CQF. We adopt 
a different perspective in our study as we don't define the hamiltonian to allow a 
comparison with IBM. 

In Sect. ([2]) we repeat a few basic definitions and we study the eigenvalues and 
eigenstates along the phase transition, in Sect. ([3]) we calculate matrix elements of the 
relevant Casimir operators and we write the potential energy functional associated to 
the phase transition, while in Sect. ^ we study the consequences of the introduction 
of higher order Casimir operators into the hamiltonian. 



2. — so(2) transition 

The book of A.Frank and P. van Isacker [B] provides a thorough pedagogical discussion 
of the u(2) dynamical symmetry and in particular of the two symmetry limits in the 
solution of the s-t-boson hamiltonian. The limits correspond to the subalgebras u{l) 
and so(2) and are associated with an anharmonic oscillator (or spherical phase) and 
with a rotor (or deformed phase) respectively. One can think of this algebra as the 
simplest possible algebra with applications to a complex system, made up of different 
interacting scalar bosons. We remind here only the essential formulas and we refer 
the reader to Ref. for further details. Notice that there are several differences in 
the definitions with respect to Ref. [7j that make the comparison a little laborious. 
We will adhere to the conventions in [5] . 

The basic building blocks are creation and annihilation operators for the two 
types of scalar bosons, called s and t. Their commutation relations are 

[s,s^^[t,t^]^l, (I) 

all the others being zero. To distinguish the two kinds of bosons one can arbitrarily 
define a parity operator, P, as follows [S]: 

Astp-i = fit Pt'fp-^ = -tt . (2) 

The following four bilinear operators, 

hg — s^s sH i^s fit = t^t (3) 

close under commutation according to the rules that define the u(2) Lie algebra. Two 
of them give directly the number operators for the boson of species s and t, while their 
sum gives the total number operator, i.e.: 

N = sKs + tH. (4) 

The bilinear operators in ^ can also be rearranged into several physically 
meaningful operators, such as, for example, the components of a angular momentum 
(see [6]). In what follows we shall only need the square of the third component of this 
angular momentum, namely: 

j2 = i(iV - s'lsht - thhs + 2sh''st) . (5) 

One can then define the most general one- and two-body hamiltonian in the s-t space 
as 

H = Eo + €nt + aril + pjj . (6) 
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In general, apart from the constant term, the hnear and quadratic terms in fit satisfy 
the m(2)Du(1) chain, while the last term satisfies the u(2)d so(2) chain. These two 
dynamical symmetry chains are the only possible chains within u{2). It is appropriate 
to rescale each of the terms in the above equation by the correct power of iV to ensure 
that the critical point will remain independent of N in the large N limit (at the leading 
order). This means that the one-body term (in this case ht) must be divided by N 
and each two-body term by iV^ as follows: 

We begin our analysis by considering a generic hamiltonian as a linear combination 
of terms with definite symmetry 

H = 6-ff„(i) + (1 - i)Hso(2) , (8) 

with Hu{i) ~ nt/N and Hso{2) = Jll^"^- The first step consists in studying the 
behaviour of the spectrum as a function of ^ for a given boson number. We do this 
here diagonalizing exactly the hamiltonian ([H]) in the u(l) basis. Notice that the only 
non-diagonal matrix elements in this basis are the ones of J\ , that are given explicitly 
in Ref. [Bj. For the sake of simplicity the spherical part has been restricted to the 
purely linear term, in order to obtain a harmonic oscillator when ^ = 1. In other 
words, Eq. ([S]) has been obtained from Eq. ([7]) setting a = 0, e = f and /3 = 1 — ^. 
The results of the numerical diagonalization are given in Fig. [1] for = 4 bosons 
and Fig. [2] for iV = 10 and N = 20. The first figure allows a direct comparison with 
the figures contained in Ref. [6j. One, of course, recovers the exact analytic solutions 
when ^ takes the two values at the limits of the interval [0, 1], corresponding to the 
two dynamical symmetry chains u{2) D so(2) and u(2) Z) At these points the 

quantum numbers coming from the two dynamical symmetries discussed above are 
strictly valid and are given on the two sides of the figure: /i labels pure so(2) states 
and Ut labels pure u{l) states. 

The energy level pattern dependence on the transition parameter for boson 
numbers 10 and 20 is displayed in Fig. [21 

This allows to compare the density of states of the two regimes (quadratic on 
the left and linear on the right side respectively) and to see that the additional 
degeneration of the so(2) shifts to the right for excited states. Notice that these 
pictures are very similar to the corresponding diagrams (for the U(3) case) of Ref. [5]. 

The diagonalization of the hamiltonian matrix given in Eq. ([H]) provides also the 
eigenvectors. We have plotted in fig. Othe squared amplitudes of the five eigenvectors 
of the case iV = 4 with respect to the components of the u(l)-basis. At the limit 
^ = 1 one clearly sees that these eigenvectors are pure in the u(l) basis, but, moving 
far from the limit, they acquire components from other states (with same parity). 
For example the ground state remains quite pure above the value ^ 0.42 that 
corresponds to the critical point and it becomes completely mixed only in the vicinity 
of the other dynamical symmetry limit. Similar considerations can be drawn also for 
the excited states, although with the due differences: for example a larger mixing is 
seen to take place for smaller values of ^ in the third and fifth states. In the present 
case the second and fourth state, being the only ones with 'negative' parity mix in a 
well-ordered fashion, just reversing the role of the highest and lowest components. 

If one repeats the calculation of the squared amplitudes for larger boson numbers 
iV, it is seen that, the critical point approaches the limit ^ = 1/2, and the eigenstates 
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Figure 1. Eigenstates of the transitional hamiltonian as a function of ^, for 
JV = 4. The two dynamical symmetries (black arrows) are recovered at the edges, 
where the labeling is also given according to ^ . The red arrow marks the critical 
point (see section 3). 




Figure 2. Same as Fig. [T] but for N=10 and N=20 bosons. The arrows mark 
the occurrence of the critical points, according to Eg. (1141 with = 0. 



keep as a largest component the one that would correspond to the spherical phase 
(u(l) limit) above that point. This fact was already noticed by lachello and several 
other authors for different cases. As an example we have plotted in fig. [Hthe same 
as the lowest panel of fig. [31 but for N = 7. In this case the critical point sits at 
= 6/13, but the eigenvectors keep as a leading component the rij = state of the 
limiting C = 1 spherical case down to about ^ ^ 0.35. This fact has an important 
consequence: although the spectrum (the energy levels) displays large modifications 
already around the critical point, the eigenvector have a little resilience to abandon 
the spherical configuration. This means that, proceeding from the spherical to the 
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Figure 3. Components, given as probabilities in the interval [0,1], of the 
transitional eigenstates as a function of ^ in the n(l) basis for the TV = 4 case. 
The pure u{l) eigenstates are indicated with full, long-dashed, short-dashed, dot- 
dashed and dotted lines for nt = 0, • • ■ , 4 respectively. Of course only states with 
the same parity appear as components of a given state. 



deformed case, all those spectroscopic properties, that depend heavily on the details of 
the wave functions (for example electromagnetic transitions, etc.), will tend to remain 
fairly similar to their 'spherical' counterpart, at least in a small interval around the 
critical point. 

3. Coherent states calculations 

The connection with a geometrical picture is obtained by resorting to the formalism 
of vector coherent states. They can be defined, in the present work, as 

\z,N)^^ ^ -^{zis^ + Z2t^)^ \0) , (9) 

where zi and Z2 are complex numbers. These states are normalized VA'^. Following a 
standard method, based on the connection of the commutator with differentiation, to 
calculate matrix elements of the operators that are relevant to the present discussion 
we obtain: 

{z,N\nt\z,N)^N (10) 

\zi\^ + 1^2 r 
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Figure 4. Components, given as mixing probabilities, of the lowest eigenstate 
as a function of ^ in the u{l) basis for the N = 7 case. Numbers from 1 to 7 
correspond to the quantum number nt of the limiting ^ = 1 case. 



, ..... \Z2\' 



{z,N\ni\z,N)^N +NiN~l) (11) 

, -9 , . N N(N - 1) (z*Z2- ziz* f , , 

{z,N\Jl\z,N) = - + ^^ (12) 



Among all the possible choices of z-y and Z2, several would cancel the two-body 
term, as it occurs for instance when z\ and zi are taken as real numbers. One choice 
that allows us to retain this term is to take z\ = \ and z^ = iy purely imaginary. 
Of course other choices are possible, provided that the matrix element don't become 
zero. Thus the energy functional associated with the hamiltonian ([7|) is given by: 



- + " + iv - + + ^) ' 

where we have further introduced a;'^ = ?;'^/(l + y^). We notice that the combination 
of and can be interpreted as a Landau potential, that leads to a second order 
phase transition. The critical point is obtained when the potential is purely quartic 
or, equivalently, by requiring the coefficient of the quadratic term to be zero. One 
obtains the following expression for the critical point: 
N - 1-a 

= (14) 

that has been used in the case a = to mark all the red arrows in the preceding 
figures. We point out here that, in the large boson limit (iV oo), the critical point 
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tends to 1/2. The presence of the quadratic Casimir of u{l) in the hamiltonian can 
modify the critical point and the whole phase transition considerably. 

4. Higher order Casimir of 

In section 4 we have considered only the linear Casimir operator of in the 

hamiltonian entering Eq. that is a = 0. This was partly due to a need for 

simplicity, but also because the fif term has a z-space realization that incorporates a z'* 
term. Conventional wisdom on quantum phase transitions and a large part of current 
literature on this topic associate (as we also have implied in the preceding sections) a 
phase with a dynamical symmetry chain. This is most often the case, because different 
symmetry chains contribute with different Casimir operators to the final energy 
functional. There might be cases in which, within a given symmetry chain, one could 
opt for including higher order terms, that, under certain choices for the parameters, 
can generate different phases, and perhaps give the same overall behaviour that was 
obtained by mixing up different subalgebra chains. Albeit extremely simple, the 
algebra discussed here allows us to discuss this general statement with a crystal-clear 
example: we take only the subalgebra chain passing through the unitary subalgebra, 
namely u{2) D u{l), with a combination of both the linear and the quadratic Casimir 
operators (therefore there is no term depending on contrarily to the previous 
section, or in other words f3 = 0) that reads: 

Hi = e^+ a% (15) 

with spectrum Ei = erit/N + ar^jN"^ . This is identical to the dynamical symmetry 
discussed in Ref. [6], sect. 1.4, save for the fact that we have dropped the constant 
term. One can again use the equations of the preceding section to determine the 
potential energy surface to be : 

E' (x) = {e ^ a)N + aN (N - \)x'^ (16) 

that has a second order critical point when e = —a,yN. It is clear that one has 
a 'spherical' minimum when e is above that threshold and a 'deformed' minimum 
otherwise. Thus two different phases can be generated within just one dynamical 
symmetry chain, without any need for mixing up with another symmetry, but just by 
including higher order terms. By higher order we mean powers of bilinear operators 
that, upon canonical ordering, might contain two-body terms, or higher order terms. 
In this case the deformed behaviour arises because of the presence of two-body terms 
coming from , rather than because of the mixing with two-body terms contained 
into J^. This case could be called anharmonic vibrator, because of the quadratic term 
in the spectrum. The foregoing discussion allows us to conclude that a 'vibrator' is 
not necessarily associated with a 'spherical' minimum although its underlying algebra 
are still u{2) and u{l), but indeed here we find a range of parameter values that allows 
for a 'deformed' minimum. In summary: when a = or when e = one only has a 
minimum in zero; when a > one may have either a minimum in zero for e > —a, 
the critical point as discussed above and a deformed minimum for e < —a. 

5. Conclusions 



The study of the u{2) algebra has been extended from the pure dynamical symmetries 
to the transitional region between them. The whole transitional path is covered by 
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constructing a hamiltonian that allows to pass with continuity from u{l) to so{2) 
by varying a parameter. Spectra and eigenstates are calculated by diagonalizing this 
hamiltonian in the u(l) basis. The ground state critical point is found by direct 
calculation of the potential energy surface within a coherent state formalism that 
must be defined with a non-standard choice of phases. 

In spite of what one might expect from several studies on quantum phase 
transitions in higher dimensional algebras, in this low-dimensional case the critical 
point cannot be associated with the e(l) symmetry. There are mainly two reasons 
for this fact. Firstly we are dealing here with a finite number of particles in an exact 
microscopic way, rather than solving a differential equation with an infinite square well 
potential in the N ~* oo limit. Therefore we cannot expect to determine a dynamical 
symmetry at the critical point as it happens for example in five dimensions in the 
solution of the Bohr Hamiltonian in nuclear physics [lOj or in two dimensions for the 
solution of the Bes equation for collective nuclear pairing [TT] . Secondly the euclidean 
group in one dimension would have a very peculiar structure since so(l) does not 
exists and rotations in one dimension essentially reduce to translation. 

Finally, a particular case including higher order operators of u{l) has been 
discussed. This allows to show that, despite common belief, second order Casimir 
operators of trigger the occurrence of a 'deformed' phase. This fact, that we 
exemplify in a schematic ^(2) model, to our knowledge has never been highlighted 
and its counterpart in more complicated algebras, such as those used in molecular 
and nuclear physics, might lead to a reconsideration of several algebraic models and 
introduce new facets in the description of various quantum systems. 
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